Weak intermolecular interactions in gas-phase NMR 
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Gas-phase NMR spectra demonstrating the effect of weak intermolecular forces on the NMR 
shielding constants of the interacting species are reported. We analyse the interaction of the molec- 
ular hydrogen isotopomers with He, Ne, and Ar, and the interaction in the He-CCh dimer. The 
same effects are studied for all these systems in the ab initio calculations. The comparison of the 
experimental and computed shielding constants is shown to depend strongly on the treatment of the 
bulk susceptibility effects, which determine in practice the pressure dependence of the experimental 
values. Best agreement of the results is obtained when the bulk susceptibility correction in rare gas 
solvents is evaluated from the analysis of the He-rare gas interactions, and when the shielding of 
deuterium in D2-rare gas systems is considered. 



I. INTRODUCTION 

The importance of intermolecular interactions in 
physics, chemistry, and biology does not need to be 
stressed. Intermolecular potentials determine the proper- 
ties of non-ideal gases, (pure) liquids, solutions, molecu- 
lar solids, and the behavior of complex molecular ensem- 
bles encountered in biological systems. They describe the 
so-called non-bonded contributions, as well as the special 
hydrogen bonding terms, that are part of the force fields 
used in simulations of processes such as enzyme-substrate 
binding, drug-receptor interactions, etc. A few examples 
showing important applications of intermolecular poten- 
tials include the Monte Carlo and molecular dynamics 
simulations of biological systems, studies of processes in 
the earth's atmosphere, or interstellar chemistry. 

Also the NMR spectra, in particular the observed 
chemical shifts, depend not only on the molecular struc- 
ture but also on the intermolecular forces. The changes 
due to the environment are difficult to interpret theo- 
retically and make the comparison of the computed and 
observed spectra unreliable. The role of the intermolec- 
ular forces is undoubtedly the largest in the condensed 
phase, and much smaller in dilute gas-phase solutions. 
Moreover, it is particularly small if we analyse a system 
where only weak van der Waals intermolecular forces play 
a significant role. In this work, we describe gas-phase 
NMR spectra for such systems, analyse the dependence 
of the observed shielding constants on the intermolecular 
forces, and present ab initio calculations which describe 
this dependence. 

Early NMR studies in the gas phase were reviewed by 
Rummens (l|, another review has been written in 1991 
by Jameson Q. However, the role of the intermolecu- 
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lar interactions in the gas phase was almost exclusively 
interpreted on the basis of binary collision gas model 
introduced by Raynes, Buckingham, and Bernstein [3J. 
In this RBB model the change in the shielding constant 
is qualitatively described as a sum of contributions due 
to the bulk susceptibility, neighbor-molecule magnetic 
anisotropy, polar effects, and van der Waals effects. At 
present, by applying state-of-the-art methods of quan- 
tum chemistry we should be able to predict accurately 
the small changes of the shielding constants due to weak 
intermolecular forces. For the first time this should be 
possible within an ab initio approach, which is in prin- 
ciple more reliable than the standard methods used to 
describe for instance the solvent effects in liquids, such 
as various polarizable continuum models based on clas- 
sical approximations. Therefore, a study of gas-phase 
model systems has a specific advantage for the compari- 
son between experiment and theory. 

Theoretical studies of the interaction-induced changes 
in the NMR parameters are scarce, and mostly restricted 
to supermolecule calculations of the interaction-induced 
shielding constants and spin-spin coupling constants; see, 
e.g. Refs. |4Hll| for typical applications. To our knowl- 
edge only one paper [J] analysed (comparing the theory 
with the numerical results) the asymptotic long-range be- 
havior of the shielding constant and its anisotropy in a 
dimer. Most of the papers reporting ab initio calculations 
of the NMR parameters that could directly be compared 
with the gas phase NMR experiment were devoted to 
studies of atom-atom interactions |a-[l2| (this is in sharp 
contrast with the electric properties of molecular com- 
plexes for which a general long-range theory and applica- 
tions to the optical and dielectric properties of gases are 
available [13l4l5j |). There are very few ab initio studies of 
the NMR effects of weak interactions between a molecule 
and an atom or two molecules. The whole property sur- 
face has been computed for the interactions in the C2H2- 
He and C2H2-H + complexes [6j, but no Boltzmann av- 
eraging has been performed; NMR properties were also 
examined I5| for the optimized geometries of other binary 



complexes of acetylene. 

Also on the experimental side not too much has been 
done. The effects of weak intermolecular interactions on 
NMR shielding of 1 H, 2 H and 3 He in gases are small and 
buried in the much larger bulk susceptibility correction, 
therefore a detailed analysis of such systems is practically 
impossible and the RBB model has mostly been used. We 
recall here that for 3 He the effect of the weak interactions 
is particularly small and difficult to observe (see for in- 
stance the study of gas-to-liquid shifts [la]). For 21 Ne 
the precision of the NMR measurements is limited be- 
cause the magnetically active isotope has a large nuclear 
quadrupole moment; for argon the only magnetically ac- 
tive 39 Ar isotope is radioactive. On the other hand, for 
129 Xe the effects are very large. They have been observed 
in the xenon dimer and nearly quantitative agreement of 
theory with experiment was reached in state-of-the-art ab 
initio calculations |l[tlll|. Also density functional theory 
(DFT) calculations for Xe-rare gas dimers yield satisfac- 
tory agreement with experimental data, see Refs. |17l[l8j . 
Most recently, the chemical shift of Xe dissolved in liquid 
benzene was studied in the calculations combining the 
DFT methods with the classical molecular dynamics [19j . 
However, there are no similar studies of atom-molecule 
systems applying well established state-of-the art wave- 
function methods and comparing the results with known 
experimental data. 

In this paper we fill this gap and report a joint ex- 
perimental and theoretical study of the gas phase shield- 
ing constants in the mixtures of atomic and molecular 
gases. We study the effects resulting from the weak in- 
teractions between a molecule and an atom in series of 
model systems: H2-He, H2-Ne and Bi2-Ar dimers and 
their deuterium-substituted isotopomers and in He-C02. 
For the selected magnetically active nuclei — 1 H, 2 H, 3 He 
and 13 C — we observe the dependence of the NMR spec- 
trum on the density of the solvent gas, which enables 
next a comparison of the ab initio and experimental re- 
sults. In the analysis of the NMR spectra we take into 
account the bulk susceptibility correction, dependent on 
the magnetizability of the medium and on the shape of 
the NMR sample [lg, [2fJ]. In the case of the weakly 
interacting systems which we study this correction dom- 
inates in the density dependence of the spectrum and its 
proper description is essential when we extract the infor- 
mation on the role of the intermolecular interactions from 
the experimental data and compare the experimental and 
computed quantities. 

The plan of this paper is as follows. We start with 
the virial expansion of the shielding constant in terms of 
the gas density and discuss all quantities needed on the 
route from the theory to a direct comparison with the 
experiment. This is thoroughly discussed in sec. [TTJ The 
details of the computational procedures adopted in the 
ab initio calculations, fitting of the interaction potential 
and shielding surfaces, and some numerical integration 
procedures will be discussed in sec. Mil The experiment 
is described in detail in sec. IIV1 The results of the mea- 



surements and calculations are reported and compared 
in sec. [V] Finally, sec. IVII concludes our paper. 



II. SHIELDING CONSTANTS IN THE 
GAS-PHASE SOLUTIONS 

For a binary mixture of a gas A, containing the nucleus 
X whose shielding a A (X) is observed, and gas B as the 
solvent, cr A (X) can be expressed as Q: 



a A {X) = a*(X) + a? A (X)p A + v?»(X)p B (1) 

where pa and ps are the densities of A and B, respec- 
tively, and a A (X) is the shielding in the zero-density 
limit. All higher terms in Eq. ([T]), which represents a 
truncated virial expansion, can safely be neglected if the 
experimental dependence of the shielding on the density 
is linear. The coefficients a AA (X) and a AB (X) are then 
the only terms responsible for the medium effects. They 
contain the bulk susceptibility corrections, o^,^ and 
°ibulk' an d the terms directly taking account of the inter- 
molecular interactions during the binary collisions of the 
A- A and A -B molecules: a A ~ A (X) and a A ~ B (X), re- 
spectively. The shielding parameters in Eq. ([1} are tem- 
perature dependent and for this reason all the present 
measurements are performed at the constant tempera- 
ture of 300 K. Moreover, in the experiments the density 
of A, pa, is always kept very low in order to eliminate 
the solute-solute molecular interactions and Eq. ((T|) can 
be simplified to: 



a A (X) = a A (X) + o-r(X)p 



where 



A-B i 



°~1 ( X ) = CT lbulk + °~i ( x ) 



(2) 



(3) 



Fig. H] displays, as an example, the dependence of the 
helium and deuteron magnetic shielding (given with re- 
spect to the isolated systems) on the density of the rare 
gas solvent in gaseous solutions. The plots in Fig. [1] are 
linear, which proves that Eq. © is a valid approximation 
and allows the determination of the a A (X) and a AB (X) 
shielding parameters. The part of the shielding constant 
a (X) which is exclusively due to pair intermolecular 
interactions between the solute and solvent molecules is 
given by o A ~ B {X). An inspection of Eqs. © and © 
shows that a A ~ B {X) can be extracted from the experi- 
mental results once the measured shielding constant be- 
comes linear in the gas density pb, and if the bulk sus- 
ceptibility correction, cr^, ulk , is known. 

In the experiment it is not easy to measure the gas 
number density p, but rather the pressure p. Therefore, 
the following form of Eq. © was used: 



a A (X)=a A (X)+a AB (X)p 



(4) 



For an ideal gas Eqs. © and (QJ are equivalent, and the 
coefficients a AB (X) and a AB (X) are inter-related by the 




density [mol/L] 



FIG. 1: The observed density-dependent 3 He shielding of 
atomic helium and 2 H shielding of deuterium in gaseous so- 
lutions (for comparison, the best estimate of the o"ibuikPAr 
contribution is shown). 



following simple expression: 

o? p B {X)=ot B {X)/k B T. (5) 

In general Eq. ([5]) is not valid since the pressure depends 
on the gas number density in a more complicated way: 



p = k B Tp + B 2 (T)p 2 + B 3 (T)p 3 



(6) 



where B 2 {T) and B 3 (T) are the second and third thermo- 
dynamic virial coefficients, respectively. Assuming that 
the NMR active molecules are infinitely diluted in the 
bath, B 2 (T) exclusively depends on the pair interactions 
between the molecules in the bath. The third virial coeffi- 
cient B 3 (T) also depends on the non-additive three-body 
interactions in the bath. We assume that we are dealing 
with an infinitely diluted solutions. In such a case we 
can assume that the concentration of the solute is very 
small and that the contribution of the partial pressure of 
the solute to the total pressure is negligible. This means 
that the thermodynamics of the system is described by 
Eq. ([6]) with the characteristic coefficients of the sol- 
vent used in the experiment, while Eq. @ describes the 
change of the shielding constant due to binary collisions 
of the NMR active molecule with the bath molecules. It 
is worth noting at this point that the virial expansions ([6]) 
and @ follow from the theory, and that in the virial ex- 
pansions the gas number density appears as the variable 
of the power series. 

A precise evaluation of cr^ ulk , the bulk susceptibility 
correction, is particularly important in the present work, 
because for all the nuclei the total change of the shield- 
ing due to intermolecular interactions in the gas phase is 
very small. We consider first the determination of erf^^ 
terms from the available xm — molar magnetic suscepti- 
bilities of gases — and applying the standard formula for 



a infinitely long cylindrical tube parallel to the external 
magnetic field [3, H3, EH : 



4-7T 
Clbulk = r- XM 



(7) 



where xm is given in ppm cgs and aibuik in ppm mL/mol. 
The macroscopic molar magnetic susceptibility xm is 
for closed shell systems proportional to the microscopic 
molecular magnetizability (1 ppm cgs corresponding to 
16.60529 xKT 30 JT~ 2 ). Equation © may be applied 
for the cylindrical geometry of the sample, and assuming 
that the molecules of the solvent do not interact. How- 
ever, this assumption is not always true, the cylinder is 
not infinite, and in such an approach various additional 
corrections are undoubtedly necessary to get realistic val- 
ues of er^uik (see for instance Seydoux et al. [16]). We 
use a different approach to determine cfbuik> which will 
be discussed in detail in section UlI El 

Theoretical determination of <r 1 ~ (X) requires two 
steps: ab initio calculations of the interaction potential 
and interaction-induced shielding constant for the binary 
complex A — B, and the average of the latter quantity 
with the Boltzmann factor depending on the interaction 
potential. The interaction potential V is given by the 
standard expression: 



V = E 



AB 



E A -E E 



(8) 



where Eab, Ea, and Eb are the energies of the colli- 
sional dimer A — B, and of the solvent (A) and solute (B) 
molecules, respectively. The interaction-induced shield- 
ing constant Q-^ B (X) is given by: 



A-B, 



(X) = a£-»(X)-v£(X), 



(9) 



where a ~ (X) and a A (X) are the shielding constants 
of the nucleus X in the dimer A — B and in the solute 
molecule A, respectively. Finally, a Y ~ (X) appearing in 
Eq. (|3]) is defined as: 



A-B 



(X) = 



<t PO ex P (-0V(ua,ub,R)) R'dRduJAduB, 

(10) 



where uja and ujb denote the two sets of the angles speci- 
fying the orientations of the monomers A and B, R is the 
distance between the centers of mass of the monomers, 
P = (^bT) _1 , ks is the Boltzmann constant, and T is 
the temperature in Kelvin. We note that in general the 
calculation of a l ~ (X) is not an easy task. For rigid 
molecules A and B it requires a six-dimensional inte- 
gration over five angles and one distance. For systems 
considered in the present paper the integral of Eq. (JTUJ) 
reduces to a two-dimensional integral that can easily be 
evaluated. The details of the computational procedures 
adopted in ab initio calculations, fitting, and numerical 
integration will be discussed in the next section. 



III. COMPUTATIONAL APPROACH 



the function: 



A. Ab initio calculations 

In all calculations the bond lengths of the interacting 
subsystems were kept fixed at their experimental geome- 
tries. We report below the results for H2 obtained with 
the H-H distance fixed at r(HH) = 1.449 an [22J. In test 
calculations we have verified that practically the same 
results are obtained using noticeably smaller values of 
r(HH). Thus, we can compare the same set of ab initio re- 
sults with the experimental data for different isotopomers 
of the hydrogen molecule. For He-C02, following the 
previous studies of the potential energy surface [23| , we 
have used r(CO)= 2.1944 an, an experimental value de- 
duced from the microwave spectra. 

All calculations of the energies and of the shielding 
constants have been performed with the coupled clus- 
ter method restricted to single, double, and noniterative 
triple excitations, CCSD(T). The NMR shielding con- 
stants and the magnetizabilities were obtained b y apply - 
ing the coupled cluster linear response theory |24l. |25| . 
Gauge-including atomic orbitals, GIAO's [2^, [27]], were 
used in all calculations of the magnetic properties, and we 
have systematically corrected all the interaction-induced 
changes in the energies and in the shielding constants by 
eliminating the basis set superposition error, i.e. all cal- 
culations for the monomers were done in the full basis 
of the dimer. We have used the d-aug-cc-pVXZ basis 
sets [28]]; d-aug-cc-pVQZ for the smallest H 2 -He system, 
and the d-aug-cc-pVTZ basis set for the larger H 2 -Ne, 
H2~Ar, and He-C02 dimers. The calculations of the en- 
ergies and shielding constants were performed using the 
ACES II 29] program, while the magnetizabilities were 
computed using the more recent CFOUR program [30j . 



B. Interaction potentials 



For two systems we used the available fitted interaction 



potential energy surfaces: for H2-Ar taken from Ref. [31 1 
and for He-C02 taken from Ref. [23|. These potentials 
were obtained from the symmetry-adapted perturbation 
theory (SAPT) calculations (see Refs. [32|, [33| for a re- 
view of the SAPT methodology and of the accuracy of 
the SAPT potentials). These potentials were shown to 
reproduce the high-resolution infrared spectra of the H2- 
Ar [H, H| and the He-C0 2 [II] van der Waals com- 
plexes. More importantly, they also reproduce very ac- 
curately the thermodynamic (pressure) virial coefficients 
36]. For other systems the interaction potential V(R, 9) 
was calculated by the supcrmolccular method according 
to Eq. ©. 

We use spherical coordinates defined with respect to 
the center of mass of the molecule. Calculations were 
performed for several angles 9 ranging from to 180° 
and for several radial distances R. For each angle 9 radial 
dependence of the interaction potential V was fitted with 



Ve(R) = e- a ^ R (A (9) - 
C 6 (9) C 8 (6) 



A 1 (9)R + A 2 {9)R 2 ) 



(11) 



R ( > 



R 8 



where a, An, A±, A2, Cq, and Cs were adjusted to fit the 
computed points at a given angle 9. We note parentheti- 
cally that odd powers of i? _1 do not appear in the long- 
range asymptotics of Eq. (fTTj) because the H 2 and CO2 
molecules are centrosymmetric. Next, interpolation was 
used to obtain the full interaction energy surface. The 
points calculated for a given radial distance R from each 
Vg fit were interpolated with a third-order polynomial in 
9. This procedure leads to a fitted/interpolated interac- 
tion energy surface V{R, 6), which was used in further 
calculations. 



C. Shielding constants 



The same technique was applied to obtain the 
Ofat B (R,0) surface. For each angle the radial depen- 
dence of a was fitted to the following function: 



JV 



S e (R)=e-<*W£ A «W Rk - E 



fe=0 



meM 



C m {0) 

R m • 



(12) 



where all the parameters appearing on the r.h.s. of the 
expression above were adjusted to fit the computed val- 
ues. 

For the hydrogen atom in H2-He, H2-Ne, and Eb-Ar a 
modification to the procedure described above was intro- 
duced. Since the interaction energy surface for these sys- 
tems is symmetric we were allowed to use symmetrized 
cr-surface cr(R, 9) calculated as an arithmetical average 
of interaction-induced shifts for both H nuclei. This 
improved the accuracy of the further integration of the 
cr int T function with the Boltzmann factor. 



D. Final integration 



To obtain the final result one has to calculate for 
the temperature of interest the Boltzmann average of 

A-Bi 



! (i?, 



A-B 



I dR d9 d(f>R 2 sm9 exp(-f3V(R,9)) a^~ B (R,9). 
Jo Jo Jo 

(13) 

Due to the axial symmetry of the considered systems 
the integration over <p gives 27r. Integrations over R nad 
9 were done numerically with the Mathematica [37J ] 
package. First, for each angle the radial integration was 
performed. The integration range [0,oo[ was substituted 
by [i?min, -Rmax] with properly defined i? min and i? max - 



The -R m in value was chosen to ensure that V(R m i n ,0) 
was positive and large enough to make the Boltzmann 
factor close to zero. The value of i? max was chosen in 
such a way that o^-Rmax, 9) was almost zero at R max , 
independent of the angle 9. This choice leads to i? m i n = 
2 ao and i? ma x = 300 ao for H2-He, for other systems the 
required integration range is smaller and within the same 
limits. The results obtained from the radial integration 
for each angle were interpolated with third-order function 
and this function was integrated over 9. 



E. Bulk susceptibilities 

To estimate the bulk susceptibility correction (BSC) 
we first used new values of magnetizabilitics obtained 
from CCSD(T) calculations. Using the d-aug-cc-pVQZ 
basis set we obtain for CO2 at the experimental geome- 
try -22.254 ppm cgs, with the basis set error estimated 
to be smaller than 0.2 ppm cgs. This value is consis- 
tent with the results of Ref. [381 ] and confirmed by new 
CCSD(T)/d-aug-cc-pCVQZ-unc calculations for the Ne 
and Ar atoms, which give -7.601 and -20.610 ppm cgs, 
also in agreement with Ref. 38] . The values of Xm and 
the corresponding bulk susceptibility effects derived from 
Eq. are given in Table fl] These values may only be 
considered as a crude approximation to the real cr^ ulk 
quantities. First, because the geometric factor is unable 
to reproduce accurately the susceptibility corrections in 
nuclear shielding; this problem was frequently discussed 
in the literature from the early days of NMR [39|, |40| . 
Secondly, we have used a special high-pressure tube, cf. 
scc. lIVl which was not spinning and this may induce non- 
negligible unknown effects. 



TABLE I: Solvent gas (B) magnetizabilities (ppm cgs) 
bulk susceptibility corrections (ppm mL/mol) 



and 





He 


Ne 


Ar 


C0 2 


Xm" 

B b 
°"lbulk 
JtS c 

B d 
"lbulk 


-1.8915 
7.923 

8.29(12) 
8.62(12) 


-7.601 
31.839 
28.79(6) 
29.56(6) 


-20.610 
86.333 
77.36(30) 
80.26(30) 


-22.254 
93.217 



"For consistency with Eq. 10 we use ppm cgs units. 
'Calculated from Eq. J7J. 
c Total effect observed in He— B interaction. 

^Determined from the total 3 He-B interaction and the computed 
interaction-induced o^ e ~ coefficients. 

Since a precise determination of the BSC value ac- 
cording to Eq. ([7]) is impossible, we have applied our 
own experimental approach to estimate the bulk sus- 
ceptibility corrections. It is well known that molecu- 
lar interactions between the atoms of rare gases disturb 
the 3 He shielding only to very small extent [16]. More- 
over, a description of such interactions is available from 



the theoretical studies of the shielding in these gas mix- 
tures [l2j. In the present work we have measured the 
density dependent 3 He shielding in helium, neon, and 



argon gases. It gave us the a{ 



HeHe 



r HeNe 



and 



r HeAr 



coefficients of Eq. ([2J, which were used next in Eq. ([3]) 
together with the interatomic interaction coefficients, to 
obtain erf bulk as af B (X) - a^ B {X). We have used 
the values of a? e - He , a^ e ~ Ne and af e " Ar , based on 
the theoretical results of Ref. [13]: -0.328, -0.776 and - 
2.901 ppm mL/mol, respectively (another available value 
of <Ji e ~ He , derived from the full configuration interac- 
tion calculations, but with a smaller basis set, is equal 
to -0.353 ppm mL/mol [8]). In this way wc determined 
the final values of the bulk susceptibility effects in the 
present NMR experiments, shown in Table [TJ Finally, we 
note that the problems related to precise determination 
of the bulk susceptibility effects are known, they have 
been recently analysed [4l], [42[ and discrepancies of the 
order of ~10% between the computed and experimental 
data have been observed [4l| . 



IV. EXPERIMENT 

The 1 H, 2 H, 3 He and 13 C NMR chemical shifts were 
measured on a Varian INOVA 500 spectrometer at 300 K 
operated at 500.61, 76.85, 381.36 and 125.88 MHz, re- 
spectively. 2 H and 13 C spectra were acquired with a 
standard two channel Varian switchable 5 mm probe, 
while 3 He and 1 H spectra in the self reconstructed helium 
probe [431 ] . Nitromethane-d3 was used for a lock system 
when the 1 H, 3 He, and 13 C NMR measurements were car- 
ried out. The 2 H experiments required a high-band lock 
operating on the proton signal of liquid tetramethylsilane 
(TMS). For this purpose a special set of coaxial glass cap- 
illaries was prepared and the same set was also used for 
the external referencing of all the chemical shifts. The 
set of capillaries contained nitromethane-d3 in the outer 
chamber and pure liquid TMS in the inner container. The 
capillaries were placed in a special non-spinning NMR 
tube which was used for all our measurements. The tube 
was made of zirconia and equipped with a metal valve for 
gas filling at high pressure (Daedalus Innovations, USA). 

The described sample setup was complex, the zirconia 
tube with the capillaries affects the external magnetic 
field, and therefore we could not apply Eq. ([7]) in our ex- 
perimental work (see also the discussion in Ref. [16]]). We 
have bypassed the problem of bulk susceptibility correc- 
tions performing analogous measurements of 3 He shield- 
ing in 4 He, Ne and Ar as gaseous solvents using exactly 
the same setup of the sample tube with the same set of 
capillaries. This series of measurements was specially de- 
signed for precise determination of the bulk susceptibility 
effects in our experiments, according to the approach dis- 
cussed in section ITlI E[ 

An efficient high-pressure system built in our labora- 
tory permitted the NMR investigations of the hydrogen 
and helium gases for a wide range of densities. As indi- 
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FIG. 2: A high pressure system for filling the zirconia NMR 
tube with a gas up to the pressure of 300 bar. 



cated in Fig. [2] the measurements in this system can be 
carried out continuously up to the total pressure of 300 
bar. All compartments were degassed when they were 
connected to the vacuum line, then a small amount of 
the solute gas was supplied from the same vacuum line 
and finally gas solvent was added and mechanically com- 
pressed. The pressure of gaseous solution was read by 
the calibrated gauge and converted into the number of 
moles following the van der Waals equation and appropri- 
ate coefficients for real gases [44]]. Gases: H2 (Air Prod- 
uct, 99.9999%), HD (Isotec, 98% D), D 2 (Isotec, 99.96%), 
3 He (Isotec, 99.96%), 4 He (Air Product, 99.9%), Ne (Air 
Product, 99.999%), Ar (Air Product, 99.9999%) and C0 2 
(Aldrich, 99.8%) from lecture bottles were used for the 
preparation of samples without further purification. 

We have performed all the measurements of the 1 H 
and 2 H shielding for the hydrogen isotopomers, H2, HD 
and D2 , as a function of the solvent density where helium, 
neon and argon were used as the solvents. Comparing the 
1 H and 2 H NMR signals from the H2 and D2 molecules 
we found that the width at the half maximum of the deu- 
terium signal is over an order of magnitude smaller than 
the same parameter of protons in H2, e.g. Ar/1/2 = 86.9 
Hz for H 2 in helium at 60 bar while Az^/2 = 6.0 Hz for D 2 
in helium at the same pressure. In practice, this means 
that the deuterium experiments deliver much more pre- 
cise data for the analysis than can be obtained from the 
1 H NMR observations of the H2 molecule. Consequently, 
we use next the 2 H NMR experimental data for compar- 
ison of the theoretical and experimental results. 

For each discussed system, measurements have been 
performed for more than 20 different solvent gas densi- 
ties. In each case, the linear fit represents well the density 
dependence of the results, with the adjusted coefficient 



of determination larger than 0.995. The experimental 
shielding constants were corrected for the gas imperfec- 
tion, and not only derived from the relations (01 and (|5|) . 



V. RESULTS AND DISCUSSION 

We begin the discussion of the effects of intermolecular 
interactions on the shielding constants with a brief sum- 
mary of the ab initio results. Three-dimensional plots of 
<7^- B (X) for the H 2 -He, 3 He-C0 2 , and 13 C0 2 -He com- 
plexes are presented in Fig. [3l Similar plots for H2-Ne 
and H2~Ar are not reported since their R and 9 depen- 
dence is nearly the same as for H2-He. An inspection of 
Fig. El shows that a in ^~ e (H) does not show any strong 
variations on R and 9. Only at very small intermolecular 
distances a stronger dependence shows up, but at these 
geometries the interaction potential is strongly repulsive, 
so the exponential Boltzmann factor is almost zero and 
these large variations of <J in 2 ~ e (H) do not contribute to 
<7 1 2 ~ e (H). Slightly more pronounced is the geometry 
dependence of the interaction-induced shielding for both 
3 He-C0 2 and 13 C0 2 -He. 



TABLE II: Calculated ab initio values of n 
(ppm mL/mol) 



A-B 

1 



T(K) 


of 2 - We (H) 


af 2 - iVe (H) 


of 2 - Ar (H) 


150 


-0.324 


-0.230 


-4.071 


200 


-0.352 


-0.281 


-4.025 


250 


-0.381 


-0.330 


-4.085 


280 


-0.398 


-0.358 


-4.144 


300 


-0.410 


-0.377 


-4.189 


320 


-0.422 


-0.396 


-4.237 


350 


-0.440 


-0.424 


-4.314 




crf e -°' U2 (He) 


ai U2 -""(C) 




150 


-6.579 


1.3050 




200 


-6.514 


1.2894 




250 


-6.548 


1.2871 




280 


-6.590 


1.2876 




300 


-6.622 


1.2883 




320 


-6.658 


1.2889 




350 


-6.714 


1.2898 





Let us now analyse the temperature dependence of the 
(7^ s coefficients calculated from Eq. ([TB")) . The results 
for all the systems are shown in Table [TTJ An inspection 
of the Table shows that the temperature effects are too 
small to be reliably determined from the experimental 
data (thus, in what follows we shall only compare the- 
oretical results with the experiment for T=300 K). The 
dependence of the computed a x ~ on the temperature T 
is almost linear. It is interesting to note that for H 2 -He, 
H 2 -Ne, H 2 -Ar, and 3 He-C02 systems o-^~ B decreases 
with T, while for 13 C02~He the opposite is found. The 
values at the lowest temperature, 150 K, do not differ 



(a) 




5 R/bohr 



(1)) 




(c) 




FIG. 3: Geometry dependence of the interaction-induced 
shielding constant for the (a) 3 He-C0 2 , (b) 13 C0 2 -He and 
(c) fib-He complexes. 



considerably from the room temperature data, suggest- 
ing that the quantum effects will start to play a noticeable 
role at still lower temperatures. 

The temperature dependence of o~ x ~ was studied ex- 
perimentally in glass samples as described earlier [43], 
but only for the 3 He-CC>2 system (rare gases like He, 
Ne and Ar could not in practice be used as solvents in 
such experiments). Unfortunately, it was not possible 
to achieve sufficient precision to perform a quantitative 
analysis of the results. In addition, a most important fac- 
tor required to determine cr( 3 He) in 3 He-C02 as a func- 



tion of the temperature — the temperature dependence of 
CO2 bulk susceptibility — is not known. 

Before we consider the comparison of the ab initio and 
experimental results, we recall that for all the systems 
the bulk susceptibility corrections are dominant. The 
observed density dependence of 2 H and 3 He shielding in 
gaseous solutions is shown in Fig. Q] (the depicted range 
of densities corresponds to the pressure of the solvent 
rare gas increasing up to 300 bar). For comparison, we 
have shown the effect of cr^ ulk for Ar. It is clear that for 
all the nuclei the total change of the shielding in the gas 
phase is very small, determined largely by erf^ ulk , and the 
precise evaluation of the bulk susceptibility corrections is 
crucial for the present study of weak molecular interac- 
tions. The BSC effect can be neglected in the analysis 
of the experimental data in two cases — when a spherical 
NMR sample is prepared or when the sample is fast spin- 
ning at the magic angle. Unfortunately, neither of these 
methods can provide accurate results for compressed gas 
at high density. 

Our final results, obtained for 300 K, are shown in 
Table HIT1 For each system, the BSC constitutes the es- 
sential part of the measured effect, thus a minor error in 
the evaluation of cibuik clearly leads to a very significant 
error in a x ~ . In particular, the standard approxima- 
tion, cibuik = — (4tt/3) xm.) is n °t sufficiently accurate. 
The error bars shown in Table U and Table IIIII do not ac- 
count for any systematic errors in the experiment, they 
represent only the errors of the linear fits to the observed 
density dependence of the results. We have considered 
systematic errors arising from a limited precision of the 
nominal reading of the absolute frequency, and errors in 
the control of the stability of the external magnetic field; 
they limit the precision of the measured shielding con- 
stants to ±0.015 ppm. However, let us recall that in 
many cases the determination of a 1 ~ (X) required two 
NMR experiments, one for the observed A—B binary sys- 
tem and one for the 3 He— B solvent. Consequently, the 
error bars are at least doubled, to ±0.030 ppm for the 
A — B system. Moreover, our experimental setup could 
slightly disturb the magnetic field as the sample was not 
spinning during the measurements; observing repeatedly 
the same samples we noticed deviations of up to 2 Hz in 
the measured frequencies. A complete analysis of these 
systematic errors, following the discussed precision of fre- 
quency measurements and possible deviations in the gas 
density inside the NMR tube, gives ±0.50 ppm mL/mol 
as an estimate of their contribution to the error bars in 
the Gi values. This estimate does not include the tab- 
ulated errors of the linear fitting of the results, and does 
not take into account the left-over errors in the analysis 
of the bulk susceptibility effects. 

The errors in the ab initio calculations are also difficult 
to estimate. The point-wise determined shielding surface 
is presumably accurate for the smallest H2-He system, 
the correlation and basis set errors becoming larger for 
the other systems. Although the following stages — fitting 
the potential and the property surfaces, followed by the 



Boltzmann average — appear to be straightforward, ap- 
proximations made in this part of the calculation con- 
tribute significantly to the final error bars. As shown in 
Fig. [21 there are regions of the shielding surface of op- 
posite contributions to the induced shielding constant, 
and therefore the final result depends heavily on a signif- 



J 



8 



icant cancellation of positive and negative contributions, 
which in turn depends on the potential surface. Follow- 
ing various test calculations we estimate that the errors 
of the computed o^ B should not exceed 15-20% of the 
discussed above final ab initio values. 



TABLE III: Measured and calculated interaction-induced effects (ppm mL/mol) 





Measured 


Estimated 


Experimental 


Calculated 




AB 




A-B 


A-B 




<J\ 


Clbulk 


a i 


G X 


<t(D) in D 2 -He 


8.07(8) 


8.62(12) 


-0.55(20) 


-0.41 


<r(D) in D 2 -Ne 


29.83(3) 


29.56(6) 


0.27(9) 


-0.38 


ct(D) in D 2 -Ar 


76.44(32) 


80.26(30) 


-3.78(62) 


-4.19 


cr(He) in 3 He-C0 2 


84.7(24) 


93.22 


-8.5(24) 


-6.62 


<r(C) in 13 C0 2 -He 


11.09(9) 


8.62(12) 


2.47(21) 


1.29 



r 



VI. CONCLUSIONS 

In this paper, we reported the first measurement of 
the changes of the NMR shielding constants due to weak 
intermolecular interactions. It became possible due to a 
new approach for the determination of bulk susceptibil- 
ity effects, which are dominant in the studied systems. 
The interpretation of the results is related to the corre- 
sponding ab initio calculations, and we observe qualita- 
tive agreement of the ab initio values with those derived 
from the experimental data. There is a series of approx- 
imations that should be analysed to improve this agree- 
ment. In particular, it is obvious that one cannot expect 
quantitative agreement without a better description of 
the bulk susceptibility effects. We have bypassed this 
problem transferring the necessary information from one 
set of the experimental data — for pairs of rare gas atoms 
systems — to another, that is to the studied molecule- 
atom systems. Such an approach appears to yield satis- 
factory results in our case, but in general a better the- 
ory, describing accurately the bulk susceptibility correc- 



tions, is needed. For larger systems, for instance involv- 
ing molecule-molecule interactions, the experiment may 
be easier, but without a proper description of these ef- 
fects the interpretation of the results is almost impossi- 
ble. Last but not least, we note that the corresponding 
theoretical calculations are also demanding, high level of 
the ab initio theory is required to obtain a reliable de- 
scription of the small changes of the shielding constants 
due to weak intermolecular forces. For larger systems it 
may be difficult to achieve satisfactory accuracy of the re- 
sults, in particular when the effects due to different parts 
of the shielding surface partially cancel out. 
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